Crack Front Waves and the dynamics of a rapidly moving crack 
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Crack front waves are localized waves that propagate along 
the leading edge of a crack. They are generated by the in- 
teraction of a crack with a localized material inhomogeneity. 
We show that front waves are nonlinear entities that trans- 
port energy, generate surface structure and lead to localized 
velocity fluctuations. Their existence locally imparts inertia, 
which is not incorporated in current theories of fracture, to 
initially "massless" cracks. This, coupled to crack instabili- 
ties, yields both inhomogeneity and scaling behavior within 
fracture surface structure. 

PACS numbers:46.50.+a, 62.20.Mk, 68.35.Ct 



dynamics of the ID front defined by the leading edge 
of the crack in a 3D material must now be considered. 
In this Letter we will demonstrate that localized waves, 
"front waves", generated by an asperity, fundamentally 
affect a crack's dynamics. Front waves (FW) are elastic 
waves that propagate along a moving crack front. We 
will show that these waves both transport energy along 
the crack front and locally impart inertia to the initially 
"massless" crack. This leads to a sharp localization of the 
energy flux at points along the front as well as both self- 
perpetuating inhomogeneities and scaling of the fracture 
surface structure. 
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Dynamic fracture is of fundamental and practical im- 
portance. We consider the behavior of a crack interacting 
with a localized defect. We shall show that this interac- 
tion can induce fundamental changes to a crack's long- 
term dynamics. These changes imply that a necessarily 
3D view of fracture must replace the basically 2D the- 
ory that is currently used to describe fracture in ideal 
materials. 

In ideal (defect- free), brittle amorphous materials, ex- 
periments [0-|| indicate that until a crack bifurcates, its 
dynamic behavior is in excellent agreement with an equa- 
tion of motion 0,p| based on a linear elastic description 
of a moving crack in a 2D material. Balancing the energy 
flux, G, per unit length of the crack with the fracture en- 
ergy, T, defined as the energy needed to create a length 
of new fracture surface yields: 



G(v,l)~G(l)(l-v/v R ) 



(1) 



where v and vr are, respectively, the instantaneous 
crack velocity and the Rayleigh wave speed of the mate- 
rial. G(l) is dependent solely on the instantaneous crack 
length, Z, and the loading conditions. As Eq. (|l|) has no 
inertial terms, a crack tip in a 2D material should be- 
have as a massless, point-like object. In ideal materials, 
Eq. (||) was shown [|| to break down when, beyond a 
critical velocity, v c , a single crack loses stability [^-|[ to 
a state in which a crack undergoes repetitive, short-lived 
microscopic branching ("micro-branching" ) events (see 

e.g. 1). 

Let us now consider the dynamics of a crack in a "non- 
ideal" material populated by asperities (i.e. defects which 
locally perturb T). When a crack encounters an asperity, 
the system's translational invariance normal to the prop- 
agation direction (z axis in Figure |l|a) is broken. Thus, 
a crack tip can no longer be idealized as a pointlike ob- 
ject propagating within an, effectively, 2D material. The 
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FIG. 1. (a) The translational symmetry along a crack front 
(z direction) is broken when the front (propagating in the x 
direction) encounters a localized asperity, (b) The interaction 
with a single asperity (located at the origin) produces local- 
ized front waves (FW) that propagate along the front while 
generating tracks on the fracture surface. Shown is a fracture 
surface scan of a 1.5mm x 1.5mm section of the xz plane. 
(c) 3 sequential profiles in the xy plane showing FW form 
and motion. (Time increases in direction of arrow.) (d) FW 
have a unique characteristic profile. Shown are 3 different FW 
whose width was normalized by the spatial scale, W, of the 
initial asperity: W = 130/xra (□) 205/ira (o) and 445/xra (A). 
Arrows indicate the propagation direction. 

FW were first predicted as propagating velocity fluctu- 
ations confined to the fracture plane (y=0 plane in Figure 
Ilia). Ramanathan and Fisher 113], building on work by 
Willis and Movchan ||ll| , |l2[ , discovered that this new type 
of elastic wave is supported by the linearized equations 
describing the perturbed stress field of a moving crack. 
They are generated by asperities and propagate at ve- 
locities, 0.94^r < Cf< vr relative to the asperity that 
produced them |lQ| , [l3|] . Thus, the FW velocity, C||, along 



the propagating front is c\\ 



FW are stable 



for T(v) — const and decay if T(v) increases with v. FW 
were also observed numerically by Morrissey and Rice 
[ |l3| , [l4[ and, under repeated interactions with asperities, 
were shown to lead to progressive roughening of the crack 
front profile, in agreement with previous predictions |15|] 
of scalar models of fracture. Resonant effects of FW were 
anticipated in Jlq ]. 

Recent experiments JlTj in glass (where T(v) ~ const 
||) have revealed that localized waves, whose propaga- 
tion velocity corresponds to the predicted FW, are indeed 
emitted when a crack interacts with an asperity. The ob- 
served waves have a distinct out-of-plane (y) component 
which leaves traces along the fracture surface (Figure [j]b). 
In addition, after an initial decay, observed FW rapidly 
converge to nondecaying long-lived waves with a unique 
characteristic profile (Figure [yd). FW scale is determined 
by the asperity size. Their shape, however, is indepen- 
dent of initial conditions. FW retain this shape upon 
collisions, sustaining, like solitons, a relative phase shift. 

Our experiments were conducted in soda-lime glass 
plates of size 38x44x0. 3cm in the x (propagation), y 
(loading) and z (sample thickness) directions, respec- 
tively (see Figure |l|a). As in |lq| , samples were loaded 
using quasistatic, Mode I loading. The crack veloc- 
ity was measured with a resolution of 50 m/s at the 
plate surfaces, i.e. the z = and z = 0.3cm = z max 
planes. Surface amplitudes were measured using a mod- 
ified Taylor-Hobson (Surtonic 3+) scanning profllometer 
with a 0.01/im resolution normal to the fracture surface. 
Asperities were externally introduced within either the 
z = or z max planes by means of scribed lines of tri- 
angular cross-section. These lines, whose scales ranged 
between 100-1000 jam, locally decreased T. T was locally 
increased when these lines were filled with Super-glue ad- 
hesive. FW were generated by both asperity types, and 
above v c = 0A2vr (vr — 3370ra/s) by micro-branching 
events. Pulse profiles (e.g. Figure [j]b,c) were constructed 
in regions of constant mean velocity. 

As Figure [y indicates, the observed FW have a sig- 
nificant out-of-plane component. Do they, in addition, 
correspond to in-plane velocity fluctuations as predicted 
in |iq JL3| , [l4l ? To ascertain this, we measured the local 
velocity of the crack front on a plate face (z max plane) 
at locations at which FW, launched from an asperity on 
the opposite face (the z = plane), reached the measure- 
ment plane. The fracture surface amplitude is compared 
to the velocity fluctuations on this plane in Figure ga. 
In Figure gb this comparison is performed for FW gener- 
ated by micro-branching events. In both cases, velocity 
fluctuations of 20-30% of the mean velocity correspond 
precisely to the arrival of the FW, as indicated by the 
surface height measurements. Moreover, the two signals 
are entirely in phase. These strong phase correlations are 
seen whenever FW are generated by externally induced 
asperities. These correlations are not trivial; a given ve- 
locity peak could be generated by either a surface pro- 



trusion or indentation. (A protrusion on one crack face 
corresponds to an indentation on the other.) 

We can estimate the normal velocity component v y ~ 
Sy/8t, using the surface amplitude, 5y, (where St = dx/v 
for a pulse of spatial extent Sx). v y is less than 1% of 
the velocity fluctuations, v x , in the propagation direction. 
Thus, the relatively small out-of-plane surface (velocity) 
variations generated by FW correspond to large local in- 
plane fluctuations of v. 

These data further indicate that FW transport signif- 
icant amounts of energy along the front since, as Eq. (H) 
indicates, changes in v directly correspond to changes in 
G. This transport of energy, due to FW, explains how ex- 
periments are able to measure velocity fluctuations when, 
generally, the measurement plane at z = or z max is situ- 
ated relatively far from the micro-branch events, located 
within the plate's interior, that initiated the fluctuations. 
Thus, the intrinsic velocity fluctuations measured in ex- 
periments are, in essence, front waves. 
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The strong correlation in both the amplitude and 
phase of the in-plane velocity measurements and the out- 
of-plane surface amplitudes indicates that these quan- 
tities can be regarded as two components of the same 
wave, rather than as two independent entities. These 
correlations suggest that significant coupling exists be- 
tween in-plane and out-of-plane stress components at 
the crack front. As no linear coupling between in-plane 
and out-of-plane stress field components exists |12|, the 
3D nature of the FW is an indication of their non-linear 
character. A further indication of FW nonlinear ity is 
their characteristic shape. This shape is obtained for 
sufficiently large perturbations |17J whereas weak per- 
turbations generate dispersive waves that rapidly decay. 

In Figure || we present the surface structure generated 
in the immediate vicinity of a single asperity that was in- 
duced at z = . In addition to the FW propagating away 
from the asperity, a number of additional FW are typi- 
cally initiated directly ahead of the asperity. These ad- 
ditional waves indicate that the crack front immediately 



ahead of the asperity retains a "memory" of the asperity's 
existence. This suggests that the crack front experiences 
inert ial effects. The highly correlated FW within these 
wave-trains may explain the strong phase correlations 
between velocity fluctuations and surface amplitudes ev- 
ident in Figure Q. 

The existence of these additional waves is surprising 
in light of the predictions of 2D fracture mechanics. As 
Eq. ([[]) indicates, a local change in the fracture energy 
should cause an immediate corresponding change in the 
instantaneous velocity of a crack. In a quasi-2D mate- 
rial, a crack tip can not be influenced by the stress field 
fluctuations generated by its interaction with an asperity. 
This is because stress variations radiate away || at the 
shear wave speed, c s , which is considerably larger than 
v. Therefore, the moment that a crack's tip passes an 
asperity's immediate vicinity, it should instantly revert 
back to its initial velocity. The wave trains in Figure 
H indicate that this does not occur. The exponentially 
decaying amplitudes of the FW within these trains are 
generated ahead of the asperity. These FW have a well- 
defined period which, like their decay length (1.8W) and 
characteristic size (~ W) scale (2.3W) with the asperity's 
initial width, W. 



+00 



5K(z)oc / (a(z') - a(z))/(z' - zfdz 1 



(2) 



where (a(z') — a(z)) is the front's deviation in the x di- 
rection at point z along the front. The local value of 
K(z) (the "stress intensity factor") is proportional to 
G(z) 1 / 2 and thereby (by Eq. ([I])) drives the local front 
velocity. Thus, any part of the front that lags (over- 
takes) another will experience an increased (decreased) 
local stress that tends to stabilize a straight front. In 
the dynamic case, delayed potentials will introduce iner- 
tial effects (a(z) =>■ a(z,t) and K{z) => K(z,t)) thereby 
giving rise to oscillatory behavior of the local stress field 
[|l4| . Numerical evidence for a local increase of the stress 
field was provided in [[L3| , where a local increase of the 
front velocity was observed directly ahead of an asperity. 
The characteristic time scale for stress oscillations is |1J] 
the time, W/c\\, in which a front wave, travelling along 

the front with velocity c\\ = J A — v 2 , traverses an as- 
perity of size W. Thus, the front immediately ahead of 
the asperity feels the effects of the oscillating stress field 

at a spatial scale of Wv/c\\ = Wv/ wc? — v 2 . This scale 

is consistent with the observed scaling of both the decay 
length (Figure [|b) and spatial periodicity (Figure ||a and 
Figure Ep) of the structure formed ahead of an asperity. 
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FIG. 3. (a) Profile (in the z = plane) of a train of FW 
of decaying amplitude generated by a single asperity. (Asper- 
ity location is between 200/xra and 400/ira). Note: FW are 
generated ahead of the initial asperity, (inset) Topographic 
map in the xz plane showing the propagating FW train. The 
maximum amplitudes (b) of the FW generated ahead of an 
asperity exponentially decay in the propagation direction, x, 
with a characteristic decay length of 1.8VK, where W is the as- 
perity width. Values of W = 520/xra (□), 500/xra (0) 430/xra 
(V) 410/xra (o) 170/nm (star) 160/nm (open □) and 130/xra 
(open A) were used. 

What is the origin of the periodicity of the FW train? 
This nontrivial behavior of the crack front ahead of the 
asperity suggests that by breaking the system's transla- 
tional invariance along the crack front (z direction), the 
initially massless crack acquires inertia. A protrusion (in- 
dentation), breaking the front's translational invariance, 
can influence (be influenced by) other parts of the front 
via stress waves. As shown in |l^] (for a static crack 
front), the local deviation, 6K(z), of the stress field in- 
tensity from that of a straight front is proportional to: 
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FIG. 4. (a) Photograph (of width 4mm) of two paral- 
lel branch-lines, (inset) FW emitted from a branch-line. 
The energy flux (G(zj) distribution (b) is highly non-uniform 
once branch-lines are formed. G(z) was estimated from the 
fracture surface presented in (a) using measurements of mi- 
cro-branch length vs micro-branch width [|20|. Branch-lines 
are formed (c) by periodic arrays of micro-branches at many 
scales. Respective photograph heights are, from top to bot- 
tom, 2.4, 0.41, 0.12, and 0.03 mm. The period (d) roughly 
scales with the branch-line width. The predicted dependence 
of the branch-line period/width ratio (line) assumes that the 
scale in x is determined by the propagation time of a FW 
across the branch-line width. Crack propagation in (a) and 
(c) was from left to right. 

Below the micro-branching instability (v < v c ) the sur- 
face structure described in Figure || is typical of that 
generated by the interaction of a crack front with an as- 



perity. Above v C: as shown in Figure Q, micro-branching 
events |17| themselves serve as FW sources since, like 
an asperity, they serve to increase the local value of the 
fracture energy. Micro-branches in glass, however, are 
not randomly dispersed throughout the fracture surface 
but, as shown in Figure Qa, are aligned along straight 
lines in the propagation direction. Their internal struc- 
ture (Figure ||c,d) indicates rough periodicity in x [ pl[ . 
As a consequence of the increased surface area created by 
the micro-branches E8[ upon branch-line formation, the 
energy flux, G, into the front is not evenly distributed 
along the front. As Figure Hp shows, the total energy 
dissipated by a branch-line at a given location, z, along 
the front can be significantly larger than in the surround- 
ing, featureless surface. This inhomogeneous distribution 
of G, which is perpetuated for the life of a branch-line, 
indicates a nonlinear focussing of energy in the z direc- 
tion that is not inherent in current theories of fracture. 
As shown in Figure |]a, multiple "branch-lines" can co- 
exist, although they have a tendency to coalesce. These 
lines [|l8| , ^l|j2^1 are either initiated spontaneously or can 
be triggered by an asperity. 

Let us consider the behavior of a crack front imme- 
diately after a micro-branch event/asperity of width W 
occurs. At the conclusion of this event the local veloc- 
ity of the crack front will momentarily "overshoot" its 
unperturbed velocity v. For v < v c , the overshoot will 
exponentially "ring down" as in Figure g. For v > v c the 
velocity overshoot will not decay, but will instead initi- 
ate another branching event, directly ahead of the first 
micro-branch/asperity. This scenario will again repeat 
itself, thereby generating yet another branching event. In 

this picture, columns of branches spaced WvjJA — v 2 

apart, in the propagation direction, will be generated. 
This scaling behavior is shown in Figure |]c. The ini- 
tial scale, W, of the branch-lines is dynamically deter- 
mined by G, as demonstrated in [|l8|,^0[ where W is a 
roughly exponential function of v. As Figure Qd indi- 
cates, this scaling behavior is indeed observed for over 
3 orders of magnitude in W. The systematic increase 
with v of micro-branch periodicity, apparent in the fig- 



ure, is consistent with the predicted v / 



7 



v 2 behav- 



ior. Note that when v ~ v c the above picture predicts 
that an asperity should initiate a branch-line. Slowly 
decaying branch-lines initiated by an asperity are indeed 



observed in |22| near v c 

In conclusion, we have demonstrated that FW both 
transport energy along a crack front and consist of cou- 
pled in-plane and out-of-plane components. The broken 
translational symmetry along the front gives rise to lo- 
cal inertia of the front. This local inertia, when coupled 
to the micro-branching instability, provides a mechanism 
for the generation of stable, directed lines of spatially pe- 
riodic micro-branches in the propagation direction. This 
picture, which yields an explanation of both branch-line 



periodicity and scaling, may provide insight on the dy- 
namic origins of fracture surface roughness p3|| . Both 
the effective focussing, by branch-lines, of energy within 
the crack front and the spontaneous birth of local iner- 
tia within a crack front point to fundamental features of 
crack dynamics that cannot be incorporated in 2D de- 
scriptions of fracture. 
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